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Abstract. We use a dynamical systems approach to analyse the tilting spa- 
tially homogeneous Bianchi models of solvable type (e.g., types VI^ and VII;,) 
with a perfect fluid and a linear barotropic 7-law equation of state. In par- 
ticular, we study the late- time behaviour of tilted Bianchi models, with an 
emphasis on the existence of equilibrium points and their stability properties. 
We briefly discuss the tilting Bianchi type V models and the late-time asymp- 
totic behaviour of irrotational Bianchi VIIq models. We prove the important 
result that for non-inflationary Bianchi type VII^ models vacuum plane-wave 
solutions are the only future attracting equilibrium points in the Bianchi type 
YTLfi invariant set. We then investigate the dynamics close to the plane-wave 
solutions in more detail, and discover some new features that arise in the 
dynamical behaviour of Bianchi cosmologies with the inclusion of tilt. We 
point out that in a tiny open set of parameter space in the type IV model 
(the loophole) there exists closed curves which act as attracting limit cycles. 
More interestingly, in the Bianchi type VII^ models there is a bifurcation in 
which a set of equilibrium points turn into closed orbits. There is a region in 
which both sets of closed curves coexist, and it appears that for the type VIIj, 
models in this region the solution curves approach a compact surface which is 
topologically a torus. 



1. Introduction 

Spatially homogeneous (SH) models are very important in contemporary cosmol- 
ogy for two main reasons. First, these models are useful in analyzing the physical 
effects of the Universe which are affected by anisotropy in the rate of expansion. 
Second, SH Bianchi universes are relatively straightforward to analyse since the 
Einstein field equations for these cosmologies reduce to ordinary differential equa- 
tions (DEs). 

A SH cosmology is said to be tilted if the fluid velocity vector is not orthogonal 
to the group orbits, otherwise the model is said to be non-tilted Pj. Tilted universes 
are of considerable interest since they contain many features of physical relevance, 
including rotation, the conditions needed for the appearance of closed time-like 
curves and some properties of quantum cosmologies. In this paper we shall assume 
that the matter content of the universe is a perfect fluid with equation of state 
p = (7 — where 7 is constant, which includes the important cases of dust and 
radiation, 7 = 1 and 7 = |, respectively. 

Much work has been done in analyzing non- tilting SH models EJIHlll]. However, 
less is known about the tilted models. Since Bianchi models with a tilted fluid ^ 
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contain up to three additional degrees of freedom, it is evident that there is an 
increase in dynamical complexity and new features emerge with the inclusion of 
tilt. 

For the class of tilted SH models, the Einstein field equations have been written 
as an autonomous DE in a number of different ways 13 |S1 E], but due to 
the complexity of the equations a detailed analysis of the dynamics is difficult. A 
subclass of models of Bianchi type V 1111 112L and models of Bianchi type II 
|14j have been studied using dynamical systems theory. In particular, a complete 
description of the dynamics of tilted SH cosmologies of Bianchi type II was given in 
|14| . with particular emphasis on the qualitative analysis of the dynamics near the 
initial singularity and at late times. It was found that for 7 < 2, the tilt destabilizes 
the Kasner solutions leading to a Mixmaster-like initial singularity, with the tilt 
being dynamically significant, while at late times the tilt becomes dynamically 
negligible unless 7 > ^ . It was also found that the tilt does not destabilize the flat 
Friedmann model, with the result that the presence of tilt increases the likelihood 
of intermediate isotropization, supporting earlier results found in tilted Bianchi V 
models [O] . 

Recently more general tilting Bianchi type universes have been studied, with 
particular emphasis on exact solutions of the Einstein equations which are attractors 
for more general classes of solution. In ^H] a local stability analysis of a wide range 
of tilting SH universes, including Bianchi class B models and the type VIo model of 
class A, were studied to elucidate their late-time behaviours. The asymptotic late- 
time behaviour of general tilted Bianchi type VIq universes was analysed in detail in 
and this was generalized to Bianchi type VIo cosmological models containing 
two tilted 7-law perfect fluids in (17| . All of the future stable equilibrium points for 
various subclasses of tilted type VIq models, as well as for the general tilted type 
VIo models, were found (and all of the bifurcation values identified; e.g., 7 = 2/3, 
7 = 10/9, 7 = 6/5 and 7 = 4/3). 

In this paper we will generalize previous work and use a dynamical systems 
approach to analyse the tilting Bianchi type VI/i and VII/i models (and their sub- 
classes) containing a 7-law perfect fluid. We shall adopt the formalism employed 
in ^B]. These Bianchi models are of general measure in the space of all SH models 
and are sufficiently general to account for many interesting phenomena. In partic- 
ular, the Bianchi type VII/j models are the most general models that contain the 
open Friedmann-Robertson- Walker (FRW) models and are consequently of special 
interest. 

The paper is organised as follows. In the next section we shall write down 
the equations of motion using the orthonormal frame formalism and discuss the 
invariant subspaces. Various gauge choices are considered, and it is emphasised 
that the 'N-gauge' is useful for studying the stability of the plane wave solutions. 
In the next section we briefly discuss models of type V and VHq . In particular we 
show that the self-similarity breaking which occurs in the non-tilted type VIIo, is 
also present for tilted, irrotational, type VIIo models. 

We then study the late-time behaviour of various tilted Bianchi models, with an 
emphasis on the existence of equilibrium points and their stability properties. In 
particular, we study the stability properties of the vacuum plane-wave spacetimes, 
and show that for 2/3 < 7 < 2 there will always be future stable plane- wave 
solutions in the set of type IV and VII/i tilted Bianchi models. Indeed, we prove 
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a theorem concerning the future evolution in the invariant type V and type VII^ 
sets, and in the case of Bianchi type VII/i models it is shown that the only future 
attracting equilibrium points for non-inflationary fluids (7 > 2/3) are the plane- 
wave solutions. We also study the form of the asymptotic tilt velocity for these 
attractors. For inflationary type fluids (0 < 7 < 2/3) the situation is simpler; we 
show that the cosmic no-hair theorem is valid in this case. 

We then study the dynamics in the neighbourhood of the vacuum plane-wave 
solutions. In the type IV models, we discover that there is a tiny region of the pa- 
rameter space in which there exists a closed orbit, which acts as an attracting limit 
cycle. Moreover, for the type VII/i models more attracting closed orbits appear. A 
numerical analysis which is presented in a companion paper |18| . shows that there 
is an open set of parameter space such that the solution curves approach a compact 
surface which is topologically a torus in the type VII/i models. 

We also study the stability of the plane waves in Bianchi type Vlh models and 
the stability of the CoUinson-French vacuum solution in the exceptional Bianchi 
type VI„i/9 models. In the final section we summarize our results, discuss possible 
early time behaviour and briefly describe future work. Some more technical results 
are presented in the appendices. 

Finally, as stressed in 0131, an important mathematical link exists between the 
various classes in the state space hierarchy. The physical state of a cosmological 
model at an instant of time is represented by a point in the state space, which is 
flnite dimensional for SH models and inflnite dimensional otherwise. This structure 
opens the possibility that the evolution of a model in one class may be approx- 
imated, over some time interval, by a model in a more special class. Thus it is 
plausible that understanding the dynamics at one level of complexity, such as for 
example the tilting Bianchi type Vlh and VII/i models, will shed light on the possi- 
ble dynamical behaviour at a higher level, such as special classes of inhomogeneous 
cosmological models. 

2. Equations of motion 

Let us consider the Bianchi models containing a tilted 7-law perfect fluid. We 
will investigate the models having an Abelian G2 subgroup. For the Bianchi models, 
these correspond to the Lie algebras. A, where A is solvable^; i.e., types I-VII/i. 
Types VIII and IX are both semi-simple and have no Abelian G2 subgroup. For 
the Bianchi cosmologies - which admit a simply transitive symmetry group acting 
on the spatial hypcrsurfaces - we can always write the line-element as 

ds^ = -dt^ + Sabi^'^i^'', 

where u}" is a triad of one-forms obeying 

2 

and C°'fj^ depend only on time and are the structure constants of the Bianchi group 
type under consideration. The structure constants C°j^ can be split into a vector 
part a^, and a trace- free part n"-'' by [3 



A Lie algebra, g, is called solvable if its derived series, defined as go = 0, 0i = [0i— l,0i— 1], 
terminates; i.e., tliere exists a k such that gj; = 0. 
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The matrix 71°'' is symmetric, and, using the Jacobi identity, aj, = (l/2)C"j^ is in 



the kernel of ; 



The time-hke vector eg is chosen orthogonal to the group orbits and is given by Gq = 
d/dt where t is the cosmological time. It is necessary to introduce a dimensionless 
time variable, r, defined by 

where H is the Hubble scalar. Following we introduce expansion-normalised 

variables and the equations of motion can be given as an autonomous system of 
DEs (see [B]). 

The equations of motion will in general depend on the choice of gauge. We 
will choose an orthonormal frame where ei points in the direction of the vector 
Ofe — \C\g^, but we will leave the remaining frame vectors 62 and 83 defined up to 
a rotation. Explicitly, this frame rotation is given by 

(2.2) 62 ~ cos(/>e2 + sin(/ie3, §3 = — sn\(f)e.2 + cos0e3, 

where 4>is a, real function. This frame rotation defines the remaining gauge freedom 
we have in our system of equations. 

Following the notation in '19', the expansion- normalised variables are now 

13 

S+ + x/3S_ %/3I]23 
73223 S+ - V3I]-. 

"0 

N + _N23 
iV23 N ~ 7V_ 

To facilitate the remaining gauge freedom we will utilize complex variables (which 
will be written in bold typeface) : 

Nx = N-+iN23, 

(2.4) El = Ei2+iSi3, ' 
A gauge transformation is then given by 

(2.5) (Nx,Sx,Si,v) ^ (e2^*N> 



Sab — 



-2S+ 
V3Si2 

V^El3 



(2.3) 



Nat = V3 



Sx = S_ + iY.23 

V ~ V2 + iV3. 



Here, the phase 4> can be a function of time and corresponds to a rotation of frame 
according to ea. H2.2|) . From the above we see that the variables N and Sx are 
spin-2 objects and Si and v are spin-1 objects under the frame rotations. The 
complex conjugates, which we will denote with an asterisk, transform similarly: 



(2.6) 



(N;,s;,st,v* 



"St 



"S*,e 



In the following, we will explicitly leave the gauge function in the equations of 
motion and we will assume that (j) is given with respect to a frame for which (j)' — 
implies i?i = where Ra is the (expansion-normalised) local angular velocity of 
a Fermi- propagated axis with respect to the triad Sa- Thus we will replace the 
function i?i with the function 0. Explicitly, we have set i?i = —</)'. At this stage it 
is important to note that all the physical variables have to be independent of the 
gauge. Hence, the objects Nx and v are themselves not physical variables; only 
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scalars constructed from the above are gauge independent (we will define what we 
mean by scalars below). 

For the expansion-normalised variables the equations of motion are: 

= [q- 2)E+ + 3|Sip - 2|Nx P + ^ {-2vl + |vp) 



(2.7) 




= (9 


(2.8) 


K 


= (q 


(2.9) 


K 


= (q 


(2.10) 


N'x 


= {q 


(2.11) 


N' 


= (9 


(2.12) 


A' 


= (9 



x/3S? - 2Nx {tA + V3N) + ^^v' 



iNx +2^/3SxiV 



The equations for the fluid are 

n 



(2.13) 17' = — {2q-(37-2) + 27^wi + [2g(7-l)-(2-7)-75]l^'} 

(2.14) = (r + 2E+)wi-2\/3Re(Siv*)-A|v|2 + \/3Im(N;v2) 

(2.15) v' = (r-S++i(/)' + Aui-iV37Vwi)v-\/3(Sx +«Nxt;i)v* 



(2.16) V = ^^Ll!l_[(3^_4)-2(7-l)A«i-5] 
where 

^ 2 1 + (7 - 1)1/2 

S2 = E^ + |SxP + |Si|2 

S = S,bc''c^ c"ca = l, ^.'' = yc^ 

= «? + |vp, 

G+ = l + {j-l)V^, 

.2.^7) y ^ [(37 - 4) - 2(7 - l)Av^] (1 - T/^) + (2 _ ^)y2^ ^ 

1 — (7 — i)y2 

These variables are subject to the constraints 
(2.18) 1 = s2 + A2 + |Nxp + r2 

717 t;i 



(2.19) = 2S+A + 2Im(S;Nx) 



G4 



(2.20) = Si(z7V- V3A)+^S*Nx +^ 

(2.21) = + 3/i(|Nx|^ -iV^) 

The parameter 7 will be assumed to be in the interval 7 € (0, 2). 

3. The state space 

The state vector can be considered as X = Sx , Si, Nx , N, A, wi, v) modulo 
the constraints eqs. (|2.19|) . (|2.20|) and (|2.21|) . The latter is the group constraint, 
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which we assume determines the parameter h and consequently is not a free variable. 
The Bianchi identities ensure that the constraints (|2.18(I - ()2.2Q(I are first integrals; 
thus satisfying the constraints on a initial hypersurface is sufficient to ensure that 
they are satisfied at all times (see also the Appendix where the expressions for their 
time derivatives are explicitly given). 

The variable fl is determined from the Hamiltonian constraint, ea. p.f 8|l . and 
can thus be eliminated. The tuple X is 12-dimensional, but the constraints reduce 
the number of independent components to 8. In the above equations the function (jj 
carries the choice of gauge and does not have an evolution equation. Specifying this 
function determines the gauge completely. We can use this function to eliminate 
one degree of freedom. This reduces the dimension of the state space to 7 (for a 
given h); i.e., the physical state space is 7-dimensional. Again we note that all of 
the physical information can be extracted from considering real scalars. 

Let us list some relevant invariant subspaces necessary for this discussion. 

(1) T{A) The full state space of tilted solvable Bianchi models. 



(2) T{Vh) Type VI„: |Nx P - TV^ > 0. 

(3) T{VIh) Type VII,.: |Nxp -N^<0. 

(4) T{VIo) Type VIq: |Nx \^ -N^ >0, A^O. 

(5) T{VIIo) Type VIIq: |Nx p - TV^ < 0, A = 0. 

(6) T{V) Type V: |Nx | = iV = 0, A 7^ 0. 

(7) T{IV) Type IV: |Nx P ~_N^ = 0, A 7^ 0. 

(8) T{II) Type II: |Nx \^ ~ = 0, A = 0. 

(9) Bil) Type I: |Nx| = ^ = A = 0. 

(10) Ti(^)Non-exceptional, one-component tilted fluid: Si = v = 0. 

(11) B{A) Non-tilted, non-exceptional: Si = v = wi = 0. 

(12) dT{A) "Tilted" vacuum boundary: Q = 0. 

(13) dB{A) Non-tilted vacuum boundary: fl = v = vi = 0. 

Note that the Hamiltonian constraint, ensures that 



In addition we will require that the tilt velocities are not superluminal; i.e. 



The variable N may or may not be bounded. However, from ea. H2.21|l we see that 
N can only be unbounded for T{VIIq). Hence, for all models except type VIIq, 
the state space is compact. 

3.1. Equilibrium points and gauge independent quantities. 

Definition 3.1 (Scalar). We will call a quantity which is independent of the trans- 
formation, cq. (|2.2I) . a scalar. Scalars are thus gauge independent quantities. 

Examples of such scalars are easy to construct. Note that S4., N and A are all 
scalars. Furthermore, objects like NxNx, NxV^, Siv*, are also scalars. In the 
above formalism it is also necessary to provide a gauge independent definition of 
equilibrium points. 



(3.1) 



E^ + |SxP + |Si|2 + |Nxp + A2 



< 1. 



(3.2) 



Definition 3.2 (Equilibrium points). A set P is said to be (an) equilibrium point(s) 
if all scalars are constants on P as functions of r. 
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We note that if different points in P have different scalars then the points may 
represent different physical configurations. The constancy as a function of r is then 
the usual requirement for equilibrium points. 

For the non-tilted, non-exceptional models, Si = v = 0, and thus from the 
gauge dependent quantities Sx and Nx there are only four real scalars. These are 
related via the identity 

(3.3) |Nx|'|Sx|'= [Rc(NxS;)]'+ [lm(NxS;)]'. 

Hence, there are only three independent scalars constructed from Sx and Nx. In 
the non-tilted analysis, Hewitt and Wainwright [2111 used these scalars as dynamical 
variables (instead of Sx and Nx themselves) at the cost of obtaining the additional 
constraint H3.3|l . In principle, we could follow the same approach for the tilted anal- 
ysis. This choice might be appropriate for certain subspaces (like Ti(A))', however, 
in general, this would lead to too many constraints (similar to ea. (|3.3fl ) to be useful 
in practice. 

3.2. Some Choices of Gauge. There are different convenient choices of gauges 
for the different invariant subspaces. Here we will discuss three such choices which 
are useful for different purposes. 

3.2.1. "F-gauge": (p' = 0. The simplest possible choice of 4> is to let = 0. Note 
that this does not specify the gauge completely; we still have a constant U{1) 
gauge freedom left. However, the F-gauge is convenient for several reasons. First, 
the system of equations are defined everywhere on the state space. Second, the 
equations of motion reduce to an autonomous system of equations which define a 
dynamical system (with 8 variables). However, in this gauge there is still a constant 
C/(l)-transformation remaining which can, for example, be used for simplifying the 
initial data. The initial data will then, modulo this constant transformation, define 
a unique path in the state space. 

It should also be noted that in this gauge the equilibrium points of the dynamical 
system do not necessarily have X' = 0. Hence, for calculating equilibrium points 
there might be a more convenient gauge choice. 

3.2.2. "N -gauge" : Nx purely imaginary. This choice of gauge uses the gauge- 
function (j) to simplify the function Nx. We can, for example, choose the function 
Nx to be purely imaginary. To ensure that Nx is purely imaginary at all times we 
have to choose 0' = V3AI]_, where A is defined by 

(3.4) 7V = AIm(Nx). 

The equation for TV can then be replaced with an equation for A which ensures a 
closed system of equations. The advantages of this gauge is that there is no gauge 
freedom left and the equilibrium points are given by X' = 0. However, care is needed 
because the parameter A is ill-defined near Nx = 0. Close to Nx — the parameter 
A will diverge and the dynamical system will cease to be well defined. Note, however, 
that in T(VIh) this cannot happen since A^ < 1 in T{VIh)- However, in T{VIIh) 
Nx = may occur. If one wants to find the equilibrium points and investigate 
their stability then this choice of gauge is particularly useful. Equilibrium points 
in the neighbourhood of Nx =0 have to be treated separately. 
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3.2.3. "Til /\'- gauge" : Si/v purely real. These choices of gauge are related and 
use the function (j) to simphfy the function Si or v. For example, we can choose 
Ini(Si) = or Ini(v) = 0. Both of these choices can be achieved by choosing (f> to 
obey 



The remaining constant cj) transformation can then be used to choose either Si or 

V purely real (however, in general both of them cannot be purely real). For this 
gauge, all of the variables are bounded and well-behaved. However, for Ti{A) the 
gauge choice becomes degenerate so that there is still an undetermined gauge for 
this invariant subspace. We observe that this invariant subspace contains all of the 
non-tilted and the vacuum models; hence, for equilibrium points on the vacuum 
boundary, or non-tilted ones, this choice of gauge is not very well suited. However, 
in T[A)—Ti{A), this choice of gauge completely fixes the gauge and is well-behaved; 
this choice was used in the tilted type H analysis |14| . 

3.3. Equilibrium points in class B spacetimes. A number of important equi- 
librium point are given in the text below. Equilibrium points in the invariant 
Bianchi type VI/j have been studied in ref. This case contains many equilib- 

rium points and we shall only present the plane-wave equilibrium points explicitly 
here (see section l53|) . 

The following result is useful in the analysis of the evolution in the invariant 
type Yllh set, T{VIIh), and the type V set, T{V): 

Theorem 3.3 (Type V and VII^ equilibrium points). All equilibrium points of type 

V (N^ ^ N ^0, A^O) and VIIh (\H^\^-N^ <Q, A^Q) have Si = v = 0; i.e., 
all equilibrium points of type V and VIIh He in the invariant subspaces Ti{V) and 
Ti{VIIh), respectively. Furthermore, equilibrium points in the interior ofTiiVIIh) 
are either plane-wave spacetimes or FRW universes. 

Proof. The type V and Vll/i cases need to be treated separately. 

Type V: Assuming v ^ leads, after a fairly straightforward calculation, to a 
contradiction. Thus we must have Si = v = 0. 

Type VIIh'. The proof is somewhat lengthy and an outline of the proof is given 
in Appendix IB. II □ 

For 7 > 2/3 the (non- vacuum) FRW universes are unstable. Therefore we im- 
mediately have: 

Corollary 3.4. In the case of Bianchi type VIIh models the only future attracting 
equilibrium points for non-inflationary fluids (7 > 2/3) are the plane-wave solu- 
tions. 

We also note that for inflationary fluids, < 7 < 2/3, the FRW universe (SI = 1, 
F = 0) is stable in T{A). In fact, for < 7 < 6/7 (!!) there exists a monotonically 
increasing function Zi, defined by 



(3.5) 



(3.6) 



7 



(l-T/2)^(2-7) 




[2(1 



7)g + (2 - 7) + 75] Zi. 
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To see that this is a monotonicaUy increasing function, note first that \S\ < 2E 
Then, using the constraint equation H2.18|l and (1 — E) > (1 — S^)/2, we can 

write 

2(l-7)g+(2-7)+75 
(3.7) > (6 - 77)S2 + 2(1 - 7)(|Nx |' + A^) + ~ ^ ^^^ r>, 

which is strictly positive for 7 < 6/7. Thus for < 7 < 6/7 is monotonicaUy 
increasing as claimed. 

Theorem 3.5. For < 7 < 6/7, all tilted Bianchi models (with ft > 0, V < IJ of 

solvable type are asymptotically non-tilted at late times. 

Proof. Use of the monotonic function Zi . □ 

In fact, the result that these models are asymptotically non-tilted at late times is 
true for all 7 < 1; this follows from an analysis, but is not covered by this theorem. 
This theorem is, in fact, true for all ever-expanding Bianchi models, including ever- 
expanding class A models^. 

An immediate result of the above is the following corollary: 

Corollary 3.6 (Cosmic no- hair). For fl > 0, V < I, and < j < 2/3 we have 
that 

lim n = l, lim F = 0. 

T — >oo r — >oo 



Proof. Using Theorem 13 . 51 we can write 

n' [2q - (37 -2) + S] 

where i5 ^ as r — s- 00. Furthermore, for < 7 < 2/3, 

2q - (37 - 2) > (2 - 37)(1 - f7) > 0, 

where equality can only occur for 17 = 1. Assuming 17 < 1, then for every < 7 < 
2/3 there will exist a ri such that 2q — (37 — 2) + S > for all r > ti. The corollary 
now follows. □ 

This result is usually referred to as the cosmic no-hair theorem; the original 
version given by Wald was for a cosmological constant. The above corollary 
generalises the non-tilted result given in P] . We shall be primarily interested in the 
future asymptotic behaviour of models with 2/3 < 7 < 2 henceforward. 



•^For class A models {A = 0) the upper bound 6/7 can be slightly improved to 14/15 by 
modifying the monotonic function Zi to 

gl-/3y/3 7 
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4. Special subsets 

4.1. Bianchi Type V. The simplest model of class B is the type V model for 
which Nx=^ = 0, ^7^0. Hewitt and Wainwright ^Sl studied a subclass of this 
model, namely the evolution in the Ti{V) invariant subspace. 
For the type V models the constraint H2.20|l simplifies to 

(4.1) VSASi = ^v, 

which implies that Si and v are parallel. Thus for type V the Si and v gauges 
are equivalent. 

Theorem 13.31 implies that all of the equilibrium points of type V are given in 
Hewitt and Wainwright's analysis jl3| . Hence, the equilibrium points can be ex- 
tracted directly from that paper. The stability analysis, on the other hand, must 
be performed in the full state space which will change the stability properties of 
the equilibrium points. 

The equilibrium points important for the late-time behaviour in the invariant 
space Ti{V) are as follows: 

(1) M, Milne: Future stable in Ti{V) for 2/3 < 7 < 4/3. 

(2) M, "Tilted Milne": Future stable in Ti{V) for 4/3 < 7 < 2. 

(3) M±, "Extremely tilted Milne": M- is future stable in Ti{V) for 6/5 < 
7 < 2. M+ unstable for < 7 < 2. 

By considering these equilibrium points in the full type T(V) state space we 
obtain the following stability properties: 

(1) M: Future stable in T{V) for 2/3 < 7 < 4/3. 

(2) M: Unstable in T{V). 

(3) M±: M- is future stable in T{V) for 6/5 < 7 < 2. M+ unstable for 
< 7 < 2. 

As we shall see later, these results can be directly extracted from the plane wave 
analysis (even though we assume Nx 7^ in this case). The reason for this is that 
the Milne universe can be considered the isotropic limit of the vacuum plane waves 
(I]+,iV--0). 



4.2. Irrotational Bianchi type VIIq. Let us take a first look at the tilted Bianchi 
type VIIo model. For the non-tilted model it is known that the variable N diverges 
as r 00 [231 which causes the non-tilted type VHq models to experience a 
self- similarity breaking at late times [231 124) . We will show that this self-similarity 
breaking persists within the subspace TiiVIlQ). In the non-tilted case, the variable 
Nx will oscillate around the value zero. To avoid the problem of the N-gauge at 
Nx = 0, we will choose the F-gauge for which cj)' = 0. 

We will consider the invariant subspace Ti{VIIq), where v = Si = 0, and we 
will assume, with no loss of generality, that iV, vi > 0. In this invariant subspace 



A DYNAMICAL SYSTEMS APPROACH TO THE TILTED BIANCHI MODELS 



11 



there are two monotonic functions, 



-|N> 



(1 + nI]^ 



(4.2) 
where 

(4.3) 



\2{l+n) 



4 



4|SxP(7V2-|NxP)(E+ + l) 



(|N> 



|Sx|2) (7V2 + |S> 



(S++n)V(l-n2)(|S.P + Il±ilf^) 



1 + nS I 



(1 _ y2)i(2-7) 



G4 



4(37-2). 



(4.5) Nx = ^{a„e-''^ + fSoe'^), = iriaoe-'^ - Poe'^ 



Theorem 4.1. For a non-inflationary perfect fluid (2/3 < j < 2) and any initial 
value in Ti{VIIo) with V < 1 and > 0, we have that 

lim |iV| = 00. 

Proof. Use of the monotonic function Z3 (the proof is completely analogous to the 
non-tilted case ^3)- '-' 

We note that this implies that these models are not asymptotically self-similar 
at late-times. Let us thus assume that 1/N —^ and investigate what consequences 
this may have for the evolution of the tilted type VIIq model. For iV — > cx), the 
equations for Nx and Sx have the asymptotic form 

N'^ = 2\/3iVSx, 

(4.4) = -2\/3iVNx. 

These equations can be solved in quadrature. By defining (p = 2\/3NdT, we 
have the asymptotic form 

i(aoe-*'^+/3oe'^), Sx=t^, 
2 Zi 

We also note that if we require Nx 7^ for all r, then |q;o| 7^ \Po\- 

We have a similar type of oscillation in this case. This can be verified by calcu- 
lating the scalar |SxP : 

(4.6) |SxP = i (|aoP + |/3o|' - 2|ao||/3o| cos2^) . 
We also note that 

(4.7) Im(SxN*x) = i(KP-|/3oP), 

so that according to the constraint (|2.19|l . the type VIIq model (A ~ 0) has non-zero 
vi if |q:o| 7^ \Po\ in the limit N — > 00. 

Based on these initial investigations, we introduce the four scalars 

= |Nx|'-l-|Sx|', a2 = 2Im(SxN;) 

(4.8) a3 = |Nx|'-|Sx|', (T4 = 2Re(SxN*x). 

At late times, the scalars and (74 will oscillate with increasing frequency. Fol- 
lowing j23|, we incorporate this oscillation into a phase ip by defining 

(4.9) (73 — psin-0, (74 — pcos-0- 
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Furthermore, by defining M = 1/N, the equations of motion can be written 
(Q - 2)E+ - ai - i- (l + S+)psinV, 



(4.1Uj 


v' 


(4.11) 




(4.12) 


v 


(4.13) 


p' 


(4.14) 


M' 


(4.15) 




where 




(4.16) 





G+ 

2(Q + E+-l)(Ti + (S+ 
V(l - V'^) 

^ ' [(37-4) 



- 1 — (Ti) psin-0, 
2S+], 



1- (7- 1)^2 

2 [(Q + S+ - l)p+ (S+ + l)(Ti sinV-- p2 gjj^^j 

-M (^Q + 2E+ - psin^A + VSMpcos?/') , 
1 

M 



4V3 + 2Af(E+ + 1)— cosV 
P 



l(37-2) + (2-7)y2 



2 G^ 

The variables and (72 are determined from the constraints, 

.2 _ i^y 



(4.17) 



f7 = 1 - (71 - 



cr2 = 



G 



+ 



Moreover, the identity H3.3|) imphes the constraint 

(4.18) al^al + p\ 

The analysis in 23_ showed that for the non-tilted model, effectively the oscillatory 
terms in the system of equations are dynamically negligible at late times. A similar 
analysis shows that the same is true in this case. This means that at late times 
we obtain a reduced system of equations by simply dropping the terms containing 
sin^/; and cosV'- This enables us to state: 

Theorem 4.2 (Global attractor for 2/3 < 7 < 4/3). For 2/3 < 7 < 4/3, a 
non-LRS universe in TiiVIIo) with < 1 and > has 

(4.19) lim (E+, (71, p,y) = (0,0,0,0), lim 17 = 1. 

r — ^00 r — >oo 

Proof. The proof is completely analogous to that in |^ (by replacing H, with /Sft, 
where /3 is defined in ea. l|4.3|l ). or using Z4 defined in Appendix IB. 21 □ 

When the fluid is stiffer than radiation (4/3 < 7 < 2), the universe will not 
isotropise in terms of the shear. By studying the reduced system of equations 
- which is obtained by simply dropping the terms containing sine and cosine in 
eas. H4.10f) - (|4.15ll ~ we can show that a solution asymptotes to a point on the fol- 
lowing line-bifurcation: 



-i(37"4), 1(37 -2), 

3(2-7) 



4(1 - 1/2) 
3(2-7) 



[(37-4)-F2(5^„4)] 



0-2 



(4.20) 



4(1-^2)^ 
0<V < 



rV(37-4)2-y2(5^_4)2_ ^ 



3V^7(2 ~ 7) 
2(1 -T/2) ' 

3(2 - 7) 



2(1 -F2) 



[l + (7-l)^'], 



37-4 
57-4" 
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Theorem 4.3 (Global attractors for 4/3 < 7 < 2). For 4/3 < 7 < 2, a non-LRS 
universe in Ti{VIIo) with V < 1 and Q > has 

(4.21) lim E+ = -i(37-4), 

and consequently approaches the line given by eg. ^4.20^ . 

Proof. A proof is provided in Appendix IB. 21 □ 

To summarise, these results imply that for perfect fluids with 2/3 < 7 < 4/3 the 
Ti{VIIo) models isotropise in terms of the shear, while for fluids with 4/3 < 7 < 2, 
the Ti{VIIo) models do not isotropise in terms of the shear. This result is similar 
to the non-tilted case [221 • Regarding the expansion- normalised Weyl tensor we 
thus expect a similar behaviour as in the non-tilted case. In the terminology of 
the Ti{VIIq) models are, for fluids stiffer than dust (7 = 1), extremely Weyl 
dominant at late times. Note that the radiation case, 7 = 4/3, has not been treated 
here; this case requires a separate analysis 

In the full tilted Bianchi VIIq case we expect N ^ 00 and a similar behaviour 
to occur. Tilted type VIIo models is presently being studied by other authors |25| . 

5. Stability of the plane waves 

In this section we will address the question: Are the vacuum plane waves stable in 
the set of tilted Bianchi models? To answer this question we must first deflne what 
we mean by the "vacuum plane waves" . These plane- wave solutions correspond to 
an invariant set, V{A), of the dynamical system. Since this set is vacuum, $1 = 0, 
the tilt velocity decouples from the equations of motion and becomes unphysical in 
the invariant set itself. However, a perturbation away from the vacuum solutions 
generically has 7^ 0, and in this case it is meaningful to talk about tilt. Hence, 
we may have a situation where the asymptotic solutions are identical but where the 
asymptotic tilt velocities are different. 

To see how the system of equations decouples near vacuum equilibrium points, 
we linearise the full system of equations with respect to an equilibrium point with 

= 0. Let Y be a column vector corresponding to the variables given by the set 
of eqs. H2.7() - (|2.12|) . and v the vector {vi,V2,V3). Furthermore, we assume that the 
constraint equations have been used to reduce the vector Y. Then, at the linearised 
level, we can write 



(5.1) 



'5Y' 




"A 


Ay 


0" 




'SY' 


sn 
















<5v 




c 


Av 


B 




Sv 



We note that the eigenvalue equation reduces to 

(5.2) det (A - All) det (B - AI2) [Xn - A) = 0. 

Hence, near the vacuum equilibrium points, the system of equations decouples and 
we can study the two decoupled systems separately. In addition, we need to check 
the value of An to make sure that the vacuum approximation remains valid. 

Since we are dealing with an invariant set for which 17 = 0, it is convenient to 
choose the Nx -gauge, where the case A = 00 correspond to N23 — 0. We recall 
that in general N = XN23. 

By vacuum plane waves we will mean all solutions where: 
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Vacuum plane waves: VIA), (0 < 7 < 2) 

g = -2S+, 7V|3 = -I]+(1 + S+), A2 = (l + E+)2, 

(5.3) = -S+(l + E+), = Si3 = S23 = = 0. 
The parameter ft, is determined by A via 

3ftS+(l - A^) = + 

and E_)_ is bounded by 

r < 0, A2 > 1 (VII,,), 

(5.4) - 1< S+ <^ < 0, A2 = 1 (IV), 

[<sjh^ A2<1(VI,). 

Without loss of generahty we can also assume that A, E_, TV > 0. 

To investigate the various versions of the vacuum plane waves (i.e., with different 
tilt velocities) we can assume that the geometric variables have the above values 
and use these in the equations for vi, V2, V3. The paradigm is that as the universe 
approaches a plane wave, the variables E+, etc. 'freeze in' and effectively 
become constants. There are several results that support this paradigm. A stability 
analysis shows that this is true locally. Furthermore, comprehensive numerical 
analysis in the full state space suggests that this is true globally; indeed, evidence 
for other behaviour has not been seen ^H]- Once these variables are frozen, we 
obtain a reduced three-dimensional dynamical system that can be studied. These 
equations can therefore be reinterpreted as describing the tilt velocity of test matter 
in a plane- wave background. For convenience, and to distinguish them from the full 

set of equations, we define (x,y,z) = (^vi , ^^^^ , ^^^^ ^ ■ The reduced equations 

for {x, y, z) are then: 

x' = (T + 2Y.+)x - A{y^ + z^) - 2V3N23yz, 
y' = {T-j:+ + Ax)y-a{l~X)z, 

(5.5) z' = {T-j:++Ax)z-a{l + \)y, 

where a = \/3(S_ — XN23) = (1 — x)\/3T,-, since A^23 = S_ for the plane-wave 
solutions. The function T is given by ea. (|2.17|l with 

=x^ +y^ + V^S = I]+(-2a;2 + y^ + z^) + 2V3E_2/z. 

The requirement < < 1, implies 

< x'^ + y'^ + z^ < I. 

Hence, the state space is the solid unit ball, D'^. 

To get some sense of the resulting dynamics, we note that 

(5.6) R' ^ (T + Ax)R, = {X + l)y^ + {X - l)z^ . 

Hence, R,x ^constant describes ellipses for type Vllh (A^ > l)i straight lines for 
type IV (A^ — 1), and hyperbolae for type VI,i (A^ < 1). (Note that since these 
equations are not U{1) invariant, there is no question that the resulting dynamics 
is a gauge effect.) 

Let us study the reduced dynamical system H5.5|l . The key observation is that 
from this reduced system, and from equation ea. H2.13|l for Q, we can say something 
about the stability of the plane waves. The strategy is as follows: We first find the 
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asymptotic behaviour of the system (|5.5|l . If, in addition, the variable fl is stable 
with respect to this asymptotic behaviour, then the vacuum plane waves will be 
stable in the full set of equations. Here, we will present heuristic arguments for 
some of the features of these model. In a companion paper |18| these features are 
more thoroughly investigated. 

The primary aim is the investigation of the type VII;i plane waves; however, it 
is very illustrative (and also convenient) to consider the type IV plane waves first. 
We will also consider the type VI^ plane-waves. 

5.1. Type IV plane waves. This case is defined by A = 1, which makes R de- 
generate. The invariant surface R = y = 0, splits the solid ball, D^, into two equal 
halves which are related via the symmetry {x,y,z) — *■ {x, —y, —z). It therefore 
suffices to consider the half in which ?; > 0. The equilibrium points are: 

(1) C{IV): Non-tilted, x^y^z^O. 

(2) C{IV): Intermediately tilted, x — 2(^Tri^(i^s+) ' 2/ = = 0, ^_^2S+ < 7 < 2. 

(3) C+ilV): Extremely tilted, x = 1, j/ = z = 0. 

(4) C-ilV): Extremely tilted, x = -1, y ^ z = 0. 

(5) ^±{IV): Intermediately tilted, x — — (i^s+')(3-^7) ' 



_ I \/(l-2S+)(4~37)(37~4-S+) _^ 4+S 
(l + S+)(3-27) 



y = Q, l^<^<„iin(|,^). 



(6) £±{IV): Extremely tilted, x = ^g^, z = ±y^lL^^I^^±l ^ y = Q 



!<!]+< -i. 



For various different values of and 7 these equilibrium points serve as attractors, 
saddle points or repellors. As we are mostly interested in the future asymptotic 
behaviour, we will only be concerned with the future attractors. In addition to the 
eigenvalues arising from the above reduced dynamical system we have to make sure 
that the eigenvalue from the f2-equation is negative as well. This will ensure that 
f2 ^ at late times so that our approximation is valid. Furthermore, one of the 
eigenvalues for jF± {IV) and £± {IV) is zero (namely, the one in the y-direction) so 
we have to go to next order to determine their stability. This is just the reflection 
of the fact that !F±{IV) and £±{IV) can only be an attractor in one of the half- 
spheres, and not in the other. 

The equilibrium points are attractors for the region y > for the following range 
of parameters: 

(1) £{IV): <7< 

(2) C{IV) : Always unstable for < 7 < 2. 

(3) C+{IV) : Always unstable for < 7 < 2. 

(4) C^{IV) : < 7 < 2, -i < S+ < 0. 

(5) T±{IV) : T-{IV) stable^ for max (||, ^^i^) < 7 < min (70, 

(6) E±{IV) : £-{IV) stable for -i < S+ < -i, < 7 < 2. 

Here, we have defined 70 for any given value of S+ as F(7o,S+) = 0, where 
F(7, E_|_) is defined in ea. (|C.3p in the Appendix. Note that there is a region where 
there are two co-existing future attractors; namely, the region where ^^^^^ < 7 < 

7o. Also, there appears to be a tiny region 70 < 7 < 5_|_2x:^ (from now on called 

(IV) is stable in the same region of 7 for the half J/ < 0. 



2- 
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Figure 1. The regions where the different equihbrium points are 
attractors for the type IV model with y > 0. In the left figure, all 
of the boundaries along the left edge mark the instability of the 
energy density, fJ. The right figure is a magnified region of the left 
one showing the loophole. 

the 'loophole') which does not have any stable equilibrium points (see Figure 0. 
This perhaps suggests that there is some interesting dynamical behaviour in this 
region. Numerical experimentation seems to indicate that there is a stable curve, 
C{IV), which acts as an attractor in this 'loophole'. This curve lies in the invariant 
subspace y — and has a 'mussel-like' shape. In Figure [21 we have plotted two 
solution curves for two different values of the parameters 7 and S+. From the 
figure one can clearly see the curve, C{IV), which acts as a limit cycle. Clearly it is 
not sufficient to only consider equilibrium points as attractors when investigating 
the dynamics close to a plane-wave spacetime in these SH models. These limit 
cycles are dealt with more elaborately in the work |18) . 

From the above we can conclude: For 2/3 < 7 < 2 there will always be future 
stable plane- wave spacetimes for the type IV model. Regarding the asymptotic tilt 
velocity, we have the following possibilities depending on the value of 7. 

• 2/3 < 7 < 6/5: The tilt is asymptotically zero i.e. non-tilted. 

• 6/5 < 7 < 4/3: The tilt can be either zero [C{IV)], intermediate [T{IV) 
or C{IV)], or extreme [C-{IV) or B{IV)]. 

• 4/3 < 7 < 2: The tilt is asymptotically extreme [C-{IV) or E(IV)]. 

5.2. Type VII/i plane waves. In this case > 1, and the level curves of R 
describe elliptic cylinders. We introduce elliptic polar coordinates by 

(VA + ly, \/A — Iz) = (i?cosv3, i?sin(p), 

where the evolution equation for ip is given by 



(5.7) 



= (l-a;)v/3(A2-l)E_. 
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Figure 2. The Mussel attractor of type IV: Two solution curves 
where the parameters are given inside (left), and just outside 
(right) the loophole. The curves are projected onto the zx-plane. 
Note the closed curve, C{IV), which acts as an attractor. On the 
right figure, the equilibrium point {IV) acts as a future attrac- 
tor. For both figures a; < and z < at late times. 

Since I]_ and A are constants for the plane waves, the tilt velocity will have a non- 
zero elliptic angular velocity with respect to the x-axis (in agreement with theorem 
. All of the equilibrium points must therefore lie on the x-axis: 

(1) C{h): Non-tilted, x ^ y = z = 0. 

(2) C{h): Intermediately tilted, x — 2(^T~i^(i'^s+) ' 2/ = ^ = 0, < 7 < 2. 

(3) C+{h): Extremely tilted, x = l,j/ = z = 0. 

(4) C-{h): Extremely tilted, x = ~\, y — z ~ 0. 

One might wonder what happened to the equilibrium points J-'±{IV) and £±{IV) 
as we went over to the type VII^ models. In fact, these equilibrium points have 
turned into closed periodic orbits. Hence, also in this case, it is not sufficient to 
consider only the equilibrium points; there are also exact closed periodic orbits that 
may act as future attractors. These closed periodic orbits are difficult to deal with 
in full generality; however, analytical, as well as numerical results concerning the 
existence of these orbits of the type Vllh model are presented in JS] . 

The stability of the type Vllh plane waves appears to be similar to that of the 
type IV model, except that the points f^±{IV) and £±{IV) have become closed 
orbits. These closed orbits appear to be stable in approximately the same region as 
their IV counterparts; however, there is some uncertainty regarding the boundary 
of the loophole. 

Regarding the closed curve C{IV) in the 'loophole', this seems to have turned 
into something more complicated for the type VII^ models. Naively, it might 
be anticipated that the attracting set C{IV) has turned into something which is 
topologically a torus, T. Numerical work confirm this claim 
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In any case we can conclude: For 2/3 < 7 < 2 there will always be future stable 
plane-wave spacetimes for the type Vllh model. Regarding the asymptotic tilt 
velocity, we have the following possibilities depending on the value of 7: 

• 2/3 < 7 < 6/5: The tilt is asymptotically zero [^(/i)]; i.e. non-tilted. 

• 6/5 < 7 < 4/3: The tilt can be either zero [/^(/i)], intermediate and oscil- 
latory, extreme [£-{h)], or extreme and oscillatory. 

• 4/3 < 7 < 2: The tilt is asymptotically extreme, and can either be non- 
oscillatory or oscillatory. 

5.3. Type VI/j plane waves. In this case we have < 1, and the level curves of 
R are hyperbolae. The invariant subspace R = defines two planes intersecting at 
y ^ z ^ and divides the ball into 4 pieces. Due to the symmetry {x, y, z) 
(x, — y, —z), the 4 pieces are 2 equivalent pairs. Regarding the equilibrium points, it 
therefore suffices to consider only equilibrium points in the region y > 0. All of the 
equilibrium points reside in the invariant subspace R — (we define h ^ 1/ \/\h\)'. 

(1) C{h): Non-tilted, x = y^z^O. 

(2) C{h): Intermediately tilted, x — 2(7^1^(1+^+) ' V ~ ~ 5+2^+ < 7 < 2. 

(3) C+{h): Extremely tilted, x — 1, y ~ z — 0. 

(4) C^{h): Extremely tilted, x = —1, y = z = 0. 

(5) ^±{h) : Intermediately tilted, x = - '^''^xS^Sf ' V = xVT^/V2, 
z = ±x\/l + A/\/2 where 

y^(l - 2I]+)[4 - 37 ± h{2 - 7)][37 - 4 - S+ T (1 + ^+)h] 



X = 



(1 + I]+)(3 - 27 ± h)^/l±h 



4 + E+±(l + E+)/i . U±2h 3 

< 7 < mm ' 



(6) E±{h) : Extremely tilted, x = 3^^+^+^|'J_^i , V = xVT^/V2, 



z 



±X\/1 + A/V2 where 



I (1 - 2E+)(-l - 4S+ T 2(1 + 

|3I]+±(l + S]+)/i| ' 

( l±2h \ 

— 1 < S+ < mm -, — 

\ 4±2h 3 + 

Here the situation is more complicated than before. However, from the eigenvalues 
we see that the equilibrium points J-'^{h) and £+{h) always have an unstable mode. 
Considering the eigenvalues for the remaining equilibrium points we note that there 
are no stable equilibrium points for h > —1 {h > 1). For h < —1 {h < 1), there 
are always some stable plane-wave equilibrium points except in a tiny region. This 
tiny region is the corresponding type Vlh loophole. As in the case of type IV, in 
the loophole there seems to be a closed curve, C{h), acting as an attractor. The 
curve C{h) lies in the invariant subspace y = xVl — A/-\/2, z = — xVl + A/V2, 
along with all of the other attractors. An illustration from the loophole is given in 
Figure 01 
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Figure 3. To the left a magnification of the type Vlh loophole 
for the plane waves (where = —h?/{i + K^)) is presented. The 
vertical line is h = \p2, — 1, which marks the transition of stability 
between C^{h) and £_{h). The right picture shows the mussel 
attractor for h = 0.41, 7 = 1.225. Two solution curves are plotted; 
one approaches C from the inside, while the other approaches C 
from the outside. 
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Figure 4. A plot of the regions of stability of the various equi- 
librium points. In the above plot we have set equal to 
the threshold value Note that there are no stable plane- 

waves for h > —\. The loophole is approximately where the line 
ft = — 3 — 2-\/2 Qi = \f2 — 1) terminates. The question marks 
indicate regions where no stable equilibrium points are presently 
known. 



By a careful analysis of the remaining equilibrium points and their eigenvalues, 

we obtain the following result: 

• For h < —\ and 7 > , there are always future stable plane- wave 

solutions. 
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• For h > —1 and < 7 < 2, there are no future stable vacuum plane- wave 
equilibrium points. 

We note that the critical unstable mode at h = —1 {h = 1) is the matter density 
n. Furthermore, we know that there are no other stable equilibrium points with 

= (for h 7^ —1/9). This is consequently a signal that for h > —1, the future 
asymptote is non- vacuum. 

In ^H] the non-tilted Collins perfect fluid solutions of type Vlh were investi- 
gated. This analysis showed that for h > —1, the Collins solutions were unstable 
against tilt whenever 7 > ^^'^'^^^^zji-' - . This threshold value is the line separating the 
regions marked "Collins perfect fluid" and those with a question mark in Figure 01 



2(3+V^) 

however, no stable tilted equilibrium points are presently known in this region. 



This instability indicates that for 7 > the late-time asymptote is tilted; 



5.4. The exceptional case: The Collinson-Ftench vacuum. The plane-wave- 
type analysis above can also be applied to the exceptional case of Bianchi type 
VI_i/g. The non-tilted analysis showed that the CoUinson- French vacuum (or the 
Robinson- Trautman solution) is stable in the non-tilted Bianchi cosmologies for 
7 > 10/9 |2ni- However, the tilted analysis showed that there is one unstable mode 
for 7 > 4/3 and one more for 7 > 14/9 0^]. Hence, a tilted analysis is required 
to check the stability. 

The Collinson- French solution is given by (in the N-gauge): 

(5.8) = E23 = = A = 0. 

In this case the Ui-equations also decouple and we can again treat these separately. 
For the exceptional case, h — —1/9, which leads to the exact vanishing of one 
of the constraint equations. Hence, the tilt-velocity can only have 2 independent 
components and we set V3 = 0. 

Writing {X,Y) — (wi,W2), we obtain 

(5.9) ¥' = (t+IV6x]y, 



where T is defined as usual with 



3 



3 



The above system possesses the following equilibrium points: 
(1) CT: X = r = 0, < 7 < 2. 



(9) CT.^- X - ^(37-4) y _ , V5(37-4)(3-27) 4 ^ ^ ^ 

Cy-i±. ^ - 2(3-7)'^-=^ V2(3-7) '3^^^2- 

(3) c:F2:X = 4g2^,r = 0,24_^<7<^. 

(4) £CT±: X = ±1, y = 0, < 7 < 2. 

We note that F = is an invariant subspace of the system, which splits the unit 
disc, , into two. These two halves are related via the symmetry transformation 
y — y. Hence, it suffices to analyse the system l|5.9|l for F > 0. 



4 / ^, / 3 
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Figure 5. Solution curves and phase portraits for the dynamical 
system H5.9|l . From top left to right: 7 = 1.2, 7 = 1.4, 7 = 1.46, 
7 = 1.48, 7 = 1.6, and 7 = 1.8. 



Based on the eigenvalues of the equilibrium points, the equilibrium points are 
future attractors in the range specified: 

(1) CT: 10/9 < 7 < 4/3. 

(2) Cji±: |<7<|+ v^. 

(3) CJ-2'- Never an attractor. 

(4) £CT-: < 7 < 2. 

In Figure El we have plotted some phase portraits for different values of 7. 

Including the previous analysis of the non-tilted equilibrium points we can thus 
conclude: 

• 2/3 < 7 < 10/9: Colhns type VLi/g perfect fluid solution is stable. 

• 7 = 10/9: Wainwright's type VP -^^g, 7 = 10/9, solution gT] is stable. 

• 10/9 < 7 < 2: The Collinson-French vacuum is stable. Moreover, for 
4/3 < 7 the asymptotic tilt velocity is non-zero and for | -I- ^'^^ < 7 it is 
extreme. 

6. Summary and Outlook 

In this paper we have considered tilted Bianchi models of solvable type. We 
have used the dynamical systems approach and written down the system of equa- 
tions with one unspecified gauge function. We pointed out that various different 
gauge choices are appropriate for different models and different applications. This 
formalism was then used to study the late-time behaviour of various tilted Bianchi 
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Bianchi 
Type 


Matter 


Attractor 


Comments 


I 


1 <7<2 




No tilt allowed 


II 


7 ' q 

7=f 
^<7<2 


C5(//) 
Hill) 
C{II) 
£{II) 


Non-tilted Collins-Stewart 
Hewitt's tilted type II [2H1 
Tilted type II bifurcation 
Extremely tilted 


IV 


^ < 7 < 2 
§ 


Plane waves 


Tilted/non-tilted 


V 


l<7<2 


Milne 


Tilted/non-tilted 


VI,,<-i 


3 ^ ' ^ l-3h 
2(1-'') ^ „ ^ o 

<7<2 


C{VIh) 
Plane waves 


Non-tilted Collins type VI^ 
i iltecl / non-tiitecl 


VI,,>_i 


2 ^ ^ 2(3+V-'i) 

3 ^ ' ^ 5+3v^ 

5+3v^ ^ / ^ ^ 


C{VIh) 

? 


Non-tilted Collins type Vlft, 
No attractors known 




a <^ < 10 
3 ^ ' ^ 9 

7-f 
i?<7<2 


C{VI_yg) 

W(F/_i/9) 

CoUinson- French 


Non-tilted Collins type VI_i/g 
Non-tilted Wainwright type Vl^i/g 
Tilted/non-tilted 


VIo 


l<7<| 
1 < 7 < 2 


C{VIo) 
C{Vh) 
£{Vh) 


Non-tilted Collins type VIq 
Tilted type VIq bifurcation 
Extremely tilted 


VII,, 


1 <7<2 


Plane waves 


Tilted/non-tilted 


VIIo 


l<7<2 


? 


Analysis incomplete 



Table 1 . The late-time behaviour of Bianchi models with a tilted 
7-law perfect fluid (see the text for details and references). The 
case 0<7<2/3is covered by Corollarv l3.6l 



models, with an emphasis on the existence of equilibrium points and their stability 
properties. 

Of particular interest are the plane-wave spacetimes of the class B models. In 
the analysis we have seen a new feature arise in the dynamical behaviour of Bianchi 
cosmologies; namely, an oscillating behaviour of the asymptotic tilt velocity. In the 
type IV model, this oscillation occurs in a tiny region of the parameter space in 
which the dynamical system undergoes a so-called Hopf-bifurcation, and is related 
to the fact that there are no stable equilibrium points for the corresponding range of 
parameters in the compact state space. This Hopf-bifurcation is a manifestation of 
a more general feature of the asymptotic behaviour of these models; their extreme 
sensitivity to the parameter values. An outline of the stability properties of the 
plane-waves was given and we pointed out that they are future attractors in the set 
of type IV and VII,, tilted Bianchi models. These models are studied further, both 
numerically and analytically, in a companion paper |18) . 

The plane-wave analysis in the type IV models has emphasised that studying 
the equilibrium points alone is not sufficient to determine their late-time behaviour. 
This becomes even more manifest for the type VII^ models where more attracting 
closed orbits appear. In addition, numerical analysis suggests that for the type 
VII,i analogue of the loophole, the solution curves will approach a compact surface 
which is topologically a torus ^H!- As noted earlier, these interesting dynamical 
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results are not gauge effects; however, since the interesting mathematical effects 
occur in the SI = invariant set, it is not clear what the physical implications of 
this behaviour is. 

In a companion paper the tilted type IV and VII/i models are studied in greater 
detail, both analytically and numerically 18 . As in this paper, we will be mostly 
interested in local behaviour close to equilibrium points, although we will also 
consider some global aspects. In particular, a more detailed analysis of the closed 
curves is performed. This analysis confirms the claims in this paper. 

The late-time behaviour of tilted Bianchi models is summarized in Tabled This 
Table is based on the results in this paper, and the following references: Type 
II [H; Type V [El; Type VI^ HSl; Type VP [Ej; Type VIq M- As can 
be seen from this table, there is a gap in our knowledge and it is of interest to 
analyze the late-time behaviour of Bianchi type Vl;i and type VIIq models. Both 
of these models seem to have a different type of behaviour to that of the other 
Bianchi models studied here; the type VI^ model possibly has a matter-dominated 
late-time behaviour and the type VIIq model might have, in the terminology of 
|29j . an extreme Weyl dominance at late times. Here we have seen that this latter 
behaviour is also present in the one-component tilted type VIIo subspace Ti(VIIq). 
The remaining Bianchi models (which are excluded in this work) are the two semi- 
simple models of type VIII and type IX. Both of these models need a different 
approach than the one adopted here. 

We have focused on the late-time behaviour in this work. An equally interesting 
task would be to consider the early time behaviour. In [Tij it was argued that the 
tilted type II is chaotic initially. We note that all Bianchi types, except for type I 
and V, have type II as part of the boundary. For the tilted type VIq model ^B] the 
type II boundary played a crucial role in the initial singular behaviour, and since 
the type II boundary is chaotic, it was argued that the tilted type VIq is chaotic 
as well. We might therefore wonder whether a similar phenomenon occurs in all 
the tilted Bianchi models having type II as part of the boundary. Indeed, since 
type V does not have the type II as part its boundary, we can ask whether tilted 
Bianchi type V model is the only tilted Bianchi model (apart from the trivial type 
I) which does not have chaotic behaviour initially. A combined dynamical systems 
and metric approach and a Hamiltonian analysis Ul suggest that this is indeed 
the case. Moreover, there is some numerical evidence for this in the current tilted 
Bianchi models under consideration. We shall return to this question in future 
work. 
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Appendix A. Consistency of the constraints 
We define the following functions: 

Qi = 2S+yl + 2Im(S;Nx) ' ^^'"^ 



Q = Ei(iiV- V3A)+iS*N> 



G+ 



G+ 

(A.l) G = A2 + 3/i(|Nxp-iV2). 

The constraints are satisfied if and only if all of these functions vanish. To check 
the consistency of these constraints, we assume first that they do not vanish and 
calculate their time-derivatives. After some algebra we obtain 

H' = 2qH + 2AQi, 

Q'l = 2(g-l + S+)Qi -2V3Re(EiQ*), 

Q' = (2(7-2-E++i0')Q-\/3SxQ*, 
(A.2) G" = 2(g + 2S+)G. 

This means that if all the constraints vanish initially, they will be zero at all times. 



Appendix B. Some proofs 

B.l. Outline of the proof of Theorem 13.31 We will prove Theorem 13. 31 for the 

Bianchi type VII/i models by contradiction. We assume, therefore, that there exists 
an equilibrium point with v ^ and TV^ _ p > and A ^ 0. The outline of 
the proof goes as follows. 

We first assume that Nx = 0. This leads to 

iV^O ^ Sx=0 ^ S?---^v2. 

Using the constraint equation we see that Vl ^ and (iTV — -x/SA)^ equals a negative 
real number. Thus A = 0, which is a contradiction. We can thus assume that 
Nx 7^ 0. This allows us to use the N-gauge, and we will assume this gauge is 
chosen henceforth. 

Let us now assume that Vt ^ 0. Using the evolution equations for A, Nx, Si, 
and V we obtain E23 = 0, T = 5, g = 2 + 35, Avi = -l-(5/2)5, 11+ = -l-(3/2)5, 
and E_ — viN2z. Inserting this information into the f2 equation we obtain the two 
separate cases: (i) V"^ = 1 , and (ii) 7 = 6/5. After some lengthy calculations 
both of these cases lead to contradictions. These calculations involve solving the 
remaining equations and making sure that all parameters are real and obey the 
required bounds; i.e., < O < 1, < < 1, etc. 

Therefore, we must have f2 = 0. In this case the equations of motion decouple 
and we can use the non-tilted analysis which implies that the only equilibrium points 
are the plane-wave spacetime. From ea. (|5.7|) we require ip' = 0, which implies x — 1 
(A^ > 1). Since x = wi, this implies that = + u| -f- u| > 1. However, since 
< 1, we must have — 1 and V2 = — which is a contradiction. Since 
V = implies Si = 0, Theorem 13 . 31 now follows. 
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B.2. Proof of Theorem 14. 3L The Proof consists of two steps: (i): The existence 
of a monotonic function, (ii): Showing that the monotonic function forces S+ — > 
— (37 — 4)/2 at late times. 

The first step is to find a monotonic function for the reduced system and amend 
it so that it is monotonic for the fuh system at late times. Such a function is: 

(B.l) ^4 = 7 ^A^TT— ^, w=-(37-4), 

^ ' (l + ml]+)2(i+™)' 2^ ' 

where /3 is given in ea. (|4.3|) . We amend this function by defining 

- Z4 P h(l - m^) - m{^+ + 1)(1 + 

l + /Afcosi/;' ~ 2\/3cri(l + mE_|_) 

Note that since \p/(Ji\ < 1 and m < 1, / is bounded. At sufficiently late times we 
have 

(B.2) S.i(^±±4 + Mi,, 

^ ' Z4 (l + ml]+) 

where B is a function of the state space variables and is bounded. Thus at late 
times, Z4 is monotonically increasing. However, for 4/3 < 7 < 2, this function is 
bounded for sufRciently late times, and hence, we get the bound 

/>oo 

(B.3) / (E+ + m)MT < cx), 

J To 

for a sufSciently large tq. 

The second step is to show that this implies S+ —rn. We will show this by 
contradiction. 

Define the everywhere positive function /(r) = (S+ +m)'^. Assume that / ^ 0. 
Then there must exist a. S > and an infinite sequence {t„} with lim„^oo r„ = 00 
such that 

(B.4) /(r„) > S, Vn. 

From the evolution equations we note that /'(r) is bounded; i.e. there exists a 
D > such that |/'(t)| < D for all r. We then define the family of functions gn 
given by 



(B.5) ff„(r) = 



(5-Z?|T„-r|, Tn-S/D <T <Tn + 6/D 

0, otherwise. 



We note that < gn{T) < /(t), and J^(7„(r)dr ~ S^/D. Moreover, by defining 
(we can without loss of generality assume that support (17^) n support (gj) — for 

00 

n 

we have < g{T) < /(r). Thus 

poo poo poo c2 

(B.6) / /(T)dr> / g{T)dT^Y. <?n(r)dT = -^l = oo, 
which contradicts ea. (jB.3|l . Hence, — > — m which proves the Theorem. 
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(C.l) 



Appendix C. Eigenvalues for non-exceptional types. 
Below, we have defined h—1/ \/\h\ and 

ih, h > 0, 
0, IV, 
h, h<0. 



0=< 



C.l. Plane waves, C{h) (for all h): 

Ai = 0, A2,3 = -2 [(1 + S+) ± 2iV3N23 \ , M = -4S+ - 
As = (37 - 4) + 2S+, A6.7 = (87 - 4) - E+ ± (1 + 

C.2. Plane waves, C.{h) (for all h): 



(37-2), 



Ai = 0, A: 



'2,3 



(1 + S+) ± 2i\/3iV23 



A4 



2-7 
7-1 



(1-2S+), 



A5 = - 



(1 - 2S]+)(57 - 6 + 27S+)(37 - 4 + 2S+) 
(^,_1)(9^_10 + 4S+(1-S+)) 

A6,7 = ii^|P^±(l-.r)(l + E,),. 
C.3. Plane waves, (for all /i): 



(l + i;+)±2iV3A^; 



Ai = 0, A2,3 = -2 
(C.2) A5 = 2(1-2E+), A6,7 = 1-2E+. 

C.4. Plane waves, £-(/i) (for all h): 



23 



A4 = 0, 



Ai = 0, A 



2,3 



(1 + S+) ± 2i^/ZN23 



-4(1 + 



^^_ 2(57-6 + 27E,) ^ A,7 = -1 - 4E, ± 2(1 + 



2-7 

C.5. Plane waves, .F±(7y): 



Ai = 0, A 



2,3 



(1 + S+) ± 2iV^N2i 



A4 = - 



(l-2S+)(57-6) 

3-27 



A5 = 0, AeAy = — 



(1 - 2E+)2(4 - 37)(37 - 4 - I]+)(3 - 27 + 7E+) 



(3-27)G(7,S+) 



Ae + A7 



(1-2S+)F(7,Sh 



4(3 - 27)(1772 - 4O7 + 24)G(7, E+) 



where 



G(7,I]+) = (57-6)S^-(18- 257 + 972)E+-3 + 27, 

F(7,E+) = [2(1772-407 + 24)2+^337-^12172-1527 + 66]^ 

(C.3) -9(7 - if (12I7* - 7367-'* + I66472 - 16567 + 612) 

Here is G(7, S+) < in the whole region under consideration. F(7, E+) = 0, 
defines a line firom (4/3,0) to (4/3, -1/4). 
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C.6. Plane waves, £±{IV): 



Ai = 0, A2 



A5 = 0, As 



(1 + E+) ± 2iV3N2 



(1-2I]+)(1+4Eh 
3S+ 



A4 



(1 -2S]+)(l + 2S^ 



2(1 -2S+)(27-3-7S^ 



3S+(2-7) 



C.7. Plane waves, T±{h): 

Ai = 0, A2,3 = -2 [(1 + S+) ± 2i\/37V23 

(l-2S+)[(57-6)±M2-7)] 
A4 = 



(C.4) 



3 - 27 ± /i 

A6A7 cx (1 - 2E+)2[(4 - 37)(1 + S+) + a(2 - 7)] 
X (37 - 4 - E+ - a) [3 - 27 + 7(E+ + a)] 
Ae + A7 = something really awful 



A5 = 2(1 - vi)a, 



where a = ±(1 + Y,^)h. The expressions for Ae and A7 are too nasty to write down, 
but Ae + A7 = defines a curve which ultimately rives rise to the type VI^ loophole. 



C.8. Plane waves, £±{h): 



Ai = 0, A2,3 = -2 
A5 = 2(l-wi)a, Ae 
A7 



(1 + S+) ± 2iV^N: 



23 



A4 



(l-2E+)(3 + 6S+ + 2a) 
3S+ + a ' 



(l-2S+)(l + 4I]+ + 2a) 
3E+ + a ' 
2(l-2I]+)[27-3-7(S++a)] 



(3E+ + a)(2-7) 
where a = ±(1 + 

Appendix D. Eigenvalues for the Collinson-French vacuum 

Below are listed the eigenvalues corresponding to the eq. (|5.9|) . Ai,2, and the 
eigenvalue corresponding to the matter density, Q, which is denoted A3. 



D.l. Collinson-French, CT: 



1 



1 



(97-14), A2 = (37-4), A3 = --(97 -10) 



D.2. Collinson-French, CTi±: 



_S± y/s^ - Ad 
Ai,2 — n ' A3 



where 



d = AiA 



1^2 



100(37 -4)(3- 27) 
3(3-7)(37+l) ' 



107(472 -97 + 6) 
'3(3-7)(372 -97+ 10) 

105 [(7- I)' 



s = Ai + A2 



721 



3(3-7)(37+l) 
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D.3. CoUinson-French, CJ-2- 

10(3-27) _ (97-14) [(97-14)^-6(7-1)^] 

3(7-1) ' ^ ^ 3(817-2587 + 202) 

107(127^-397 + 32) 
3(7- l)(8l72 - 2467+ 190)' 



D.4. Collinson- French, ECT±: 

2 / r-\ 2 |"97 - 14 + ^/6(7 - 1)1 1 ^ 

^^ = 3(1^^)' ^ 3(2-7) ' A3 = --(4+ 76). 
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